Cusp Kernels for Velocity-Changing Collisions 
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We introduce an analytical kernel, the "cusp" kernel, to model the effects of velocity-changing col- 
lisions on optically pumped atoms in low-pressure buffer gasses. Like the widely used Keilson-Storer 
kernel [J. Keilson and J. E. Storer, Q. Appl. Math. 10, 243 (1952)], cusp kernels are characterized 
by a single parameter and preserve a Maxwellian velocity distribution. Cusp kernels and their su- 
perpositions are more useful than KS kernels because they are more similar to real kernels inferred 
from measurements or theory, and are easier to invert to find steady-state velocity distributions. 

PACS numbers: 34.20.Cf, 32.80.Xx, 42.62.Fi, 95.75.Qr 



Collision kernels W(x, y) are often used to model 
the effects of binary, velocity- changing collisions in ex- 
periments with dilute, optically pumped atoms in low- 
pressure buffer gases [l|-|3|. Collision kernels are con- 
nected to the transport coefficients of classical transport 
theory 0, H[ . By "collision" we mean a statistical ensem- 
ble of collisions with impact parameters b less than some 
maximum value b m , above which the collisional effects 
are negligible, and including all orbital planes — or the 
partial wave equivalent. If a pumped atom has a velocity 
y along the direction of a pumping beam before a colli- 
sion, the post-collision probability to find the atom with 
a velocity between x and x + dx along the same direction 
is dxW(x,y). We express these velocities in units of the 
most-probable three-dimensional speed Urj = ^/2ksT/M 
for atoms of mass M and temperature T, where ks is 
Boltzmann's constant. One of the most important ap- 
plications of this physics is the backscattering of light 
by sodium atoms at an altitude of 90-100 km above the 
Earth in laser guidestar systems, which are extensively 
used in modern ground-based telescopes to compensate 
for atmospheric turbulence [H, 0] ■ 

For many years, models of experiments have most often 
used the Keilson-Storer (KS) collision kernel [8], which 
has the simple analytical form 



W a (x,y) = 



(1) 



The KS kernel has a single, real "memory parameter" a, 
with < a < 1 and a corresponding width b = yl — a 2 . 
The KS kernel is normalized as a probability distribution, 



W a (x, y)dx = 1, 



(2) 



and has no effect on a Maxwellian velocity distribution, 
e~ x /y/n, since 



W a (x,y)e y 'dy = e 



(3) 



Collision kernels that do not satisfy the constraint ([3]) 
produce non-Maxwellian steady-state velocity distribu- 
tions, which are not states of maximum entropy. 



KS kernels, however, are only modestly similar to real 
kernels measured by experiments or modeled from realis- 
tic interatomic potentials The most striking dif- 
ference is that real kernels display a sharp peak near the 
initial velocity y, due to weak collisions (or small-angle 
scattering), that is absent in KS kernels. More realistic 
model kernels with sharp peaks often lack the simplic- 
ity of KS kernels or violate the equilibrium requirement 
d3]) [15rll8| . Here, we discuss a "cusp" kernel C s (x,y), 
a model kernel parameterized by a single "sharpness pa- 
rameter" s, that satisfies ([2]) and ([3|), but is more similar 
to real kernels than the KS kernel, and is much easier to 
invert to find steady-state velocity distributions. 

Let p(x, t)dx be the number of atoms at time t with 
velocity between x and x + dx m a spin-polarized mode of 
a density matrix. A variant of the Boltzmann equation 
describes the evolution of p by 

d 

— p(x, t)=- (7 sd + 7vd)p(a:, t) 

+ 7vd / W{x, y)p{y, t)dy + P(x, t). (4) 

Collisions cause the spin polarization to relax at the 
velocity- independent rate j s d- For population spin 
modes, like longitudinal spin polarization, 7 sc j is a non- 
negative real number. For coherence spin modes, like 
transverse spin polarization, the Bohr frequency of the 
mode is included as an imaginary part of 7 S( j. Velocity- 
changing collisions transfer atoms with pre-collision ve- 
locities between y and y + dy to atoms with velocity x at 
the rate ^dWix^y). Velocity-selected spin polarization 
is produced by optical pumping at a rate P{x,t). For 
a monochromatic pumping laser one can approximate 
P{x,t) m pS(x — xi) where xi is the velocity of atoms 
that have been Doppler-shifted into resonance with the 
laser light, and p parameterizes the pumping rate. 

Snider [3] has shown that the KS kernel ([lj can be 
written in terms of its right and left eigenfunctions as 



w a (x, y ) = J2" n MxK(y)- 



(5) 



71 = 



2 



The right and left eigenfunctions are 

H n (x)e~* 2 



v n (x) 
and v£(y) 



H n {y) 



(6) 
(7) 



respectively, where H n {x) is a Hermite polynomial: 
H Q (x) = l,Hi(x) = 2x,H 2 (x) = Ax 2 - 2,... 
The right and left eigenfunctions are orthonormal, 
/ v%(x)v m (x)dx = 5 nm , and complete, Y, n v n( x ) v n(y) = 
S(x - y). 

In reality, there are many more "weak" (or grazing- 
incidence) collisions, corresponding to a KS memory pa- 
rameter awl, than "strong" (or head-on) collisions with 
a w 0. Simple KS kernels do not capture this dominance 
of weak collisions. A better model would be a superpo- 
sition of KS kernels weighted to have the most probable 
value of a close to 1. Such a probability density is 



P s (a) 



sa 



(8) 



where we will call the nonzero parameter s the "sharp- 
ness." For large positive sharpnesses, s » 1, the prob- 
ability density ([5]) heavily weights memory parameters 
awl. We call the kernel produced by the probability 
density © and ©, 



(9) 



C s (x,y) = I W a (x,y)P s (a)da 



E 

71 = 



s + n 



a cusp kernel. Cusp kernels are normalized, as in ^ , and 
they transform a Maxwellian distribution into itself, as 
in ([3]). Representative KS and cusp kernels are shown in 
Fig. [TJ Cusp kernels often resemble hard-sphere kernels 
2CH22J . For s > 1 , cusp kernels display a sharp peak near 
the initial velocity y, with a sharpness that increases with 
the value of the parameter s. 

Morgan and Happer [23j have shown that (JSJ can be 
written as 

C s {x, y) = ^^l e y2 Rsi-x^R^), (10) 
v 7r 

where T(s) is the Euler gamma function, x> is the greater 
of the two arguments x and y, and x < is the lesser. For in- 
teger values of the sharpness s, the right functions R s (z) 
can be written as 



Rn{z) 



H erfc„_i(z), n = 0, +1, +2, . . 
z2 H H (z), n = 0,-l,-2,.. 



(11) 



Here, erfc„(z) denotes a repeated integral of the error 
function, as discussed in Sec. 7.2 of Abramowitz and Ste- 
gun |24j |. For noninteger s, the right function R s (z) can 
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FIG. 1. (color online) A comparison of KS kernels (top) with 
various memory parameters a and cusp kernels (bottom) with 
various sharpness parameters s. In contrast to KS kernels, 
cusp kernels display a sharp peak near the initial velocity 
y = l. 



be evaluated with the power series 



^n!2-«r(i + §-§) 



(12) 



R s (z) is an entire function of both z and s. Using asymp- 
totic expressions for R s [231 ] . we find an approximation 
for C s (x, y) that gives values almost the same as those of 
(fTU| for \s\ S> 1 and for both \x\ and \y\ of order 1, 

2\nC s (x,y) ~ y 2 - x 2 - 2\x - y\V2~s + \n ^ . (13) 



The inverse of the transformation 



W a (x,y) 



~-C s (x,y)e 



is 

-s In a 



ds, (14) 



where c > is any positive number. To prove (|14[) . one 
can represent C s (x, y) with the series ©, close the path 
of integration with an "infinite semicircle" in the negative 
half of the complex s plane (for < a < 1), and use 
Cauchy's residue theorem to recover the series (|5|). 

For either KS kernels or cusp kernels and for time- 
independent pumping we can write Q as 

d f 

-p(x,t) = - / T(x,y)p(y,t)dy + P(x), 



dt' 



where the damping kernel is 



71=0 



(15) 



(16) 



3 



and the eigenvalues of the damping kernel for KS and 
cusp kernels are 



In 



7sd + 7vd(l-a") KS, 
7sd + 7vd?V (n + s) cusp. 



(17) 



Most experiments actually measure the steady-state den- 
sity matrix, the solution to (IT5|) when dp/dt — 0, which 
has the amplitude 



p(x) = j T-\x,y)P(y)dy. (18) 
Using (|16p to evaluate the resolvent T~ 1 (x,y), we find 

v n(x)Vn(v) _ 1 , lvdW(x,y) 



71=0 



In 



7o 



7oo7o 



i (19) 



where 70 = 7 s d and 700 = 7 S( j + 7vd- We see from ((18)) 
that the first term of 1(15)1. I/700, which corresponds to 
no velocity-changing collisions, "imprints" the primary 
velocity distribution P from the laser source onto the 
steady-state velocity distribution p. The second term 
gives a collisional background (or pedestal) that is spread 
out over additional velocities, as described by the "resol- 
vent kernel" W(x,y). For KS and cusp collision kernels 
the resolvent kernels are, respectively, 



k=0 



(20) 



and W{x,y) = C r (x,y), where r = 5(70/700)- (21) 

If we represent the collision kernel with a cusp kernel 
C s (x,y) of sharpness s, then the resolvent kernel ((21)) 
is simply another cusp kernel C r (x,y), defined by ((9]) 
or ((TO)) , with a diminished sharpness r = 5(70/700)- In 
contrast, if the collision kernel is represented with a KS 
kernel, the resolvent kernel is an infinite series (f2TJ)) of KS 
kernels with the sequence of modified memory param- 
eters a, a 2 , a 3 ,..., which converges rather slowly unless 
the parameter a is nearly zero. 

To better approximate a real collision kernel, one 
could superpose m > 1 cusp kernels with sharp- 
nesses si, S2, • • • , s m , and with corresponding weights 
fi,fs, ■ ■ ■ , fm that sum to unity, J2j fj = h as 



(22) 



The multi-cusp kernel ((22)) satisfies the constraints ((2]) 
and J3]). The second line of ((T7)) is then 



7n = 7sd + 7vd 



fjn _ N(n) 



n + Sj D(n) 



(23) 



We can write the denominator and numerator of ()23[) as 
the polynomials 

D(n) = (n + 8i)(n + s 2 ) ■■■(n + s m ) (24) 
and N(n) — joo(n + ri)(n + r 2 ) • • ■ (n + r m ), (25) 

where we have introduced the roots 77 of N(n). Substi- 
tuting ([2"3"f into ((T5)l and expanding D(n)/N(n) in partial 
fractions, we find that the resolvent kernel is 



W(x,y) = ^g 3 C rj {x,y). 



(26) 



Here, C r .(x,y) is a cusp kernel of sharpness fy, and the 
weight coefficients gj, which sum to unity, y7 ^ = 1, are 



.9.7 



lim 



7oo7o( Tt + r ]) D (-fj) 



(27) 
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FIG. 2. (color online) A comparison of a three-cusp fit and 
a measured collision kernel for Rb in He from Gibble and 
Gallagher (GG) [HJ f° r the initial velocities y = 0.06 (top) 
and y = 1.37 (bottom). The GG kernel is represented by 
a model form with parameters from Table I of [13 ) . normal- 
ized according to @. The parameterized model form is es- 
timated to agree with the measured kernel for Rb in He to 
within roughly 10% over the range of velocities shown [T3 ] . 
The GG kernel and three-cusp kernel are solid lines, with the 
GG kernel (black) slightly sharper than the three-cusp ker- 
nel (red). The three-cusp kernel was fit to the GG kernel 
for y = 0.06 (top), which gave the superposition parameters 
[A, /a, fs] = [0.13,0.37,0.50] and [si,s 2 ,s 3 ] = [7.8,27.2,500], 
as defined by (|22[) . The individual cusps fjC Sj (x, 0.06) are 
shown by the dashed lines in the top panel. As the bottom 
panel shows, the same fit parameters give good a gree ment 
with measurements for y = 1.37. See Fig. 10 of GG [l^] f° r a 
comparison with KS fits to measurements. 
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We see that if the damping kernel T(x, y) contains m 
cusp kernels of sharpnesses si, S2, • • • , s m , the resolvent 
T~ 1 (x,y) also contains m cusp kernels with modified 
sharpnesses ri,r 2 , ■ ■ ■ ,r m given by (f^5|) . There does not 
seem to be an analogous, simple relation for the resolvent 
of a damping kernel with multiple KS kernels. An exam- 
ple of how well a sum of cusp kernels can approximate 
an experimentally inferred kernel 12| is shown in Fig. 2. 
As the figure demonstrates, cusp kernels are a convenient 
basis to parameterize experimentally inferred kernels. 

In summary, we have introduced a convenient model 
collision kernel, the cusp kernel C s {x,y) 1 which is char- 
acterized by a single sharpness parameter s. Like KS 
kernels, cusp kernels and their superpositions are nor- 
malized, as in ©, and they have a Maxwellian velocity 
distribution for their equilibrium state, as in ([3]). Com- 
pared to KS kernels, cusp kernels and their superposi- 
tions are more similar to real kernels and can be more 
conveniently inverted to model the steady-state velocity 
distributions of optically pumped atoms in low-pressure 
buffer gases. 

The authors are grateful to Paul R. Berman, Dima 
Budker, Kurt E. Gibble, and A. M. Shalagin for helpful 
discussion. This work was supported by the Air Force 
Office of Scientific Research. 
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